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Abstract 

We analyze the space of linearized non-supersymmetric deformations around a IIA 
solution found by Cvetic, Gibbons, Lii and Pope (CGLP) in hep-th/0101096. We 
impose boundary conditions aimed at singling out among those perturbations those de- 
scribing the backreaction of anti-D2 branes on the CGLP background. The correspond- 
ing supergravity solution is a would-be dual to a metastable super symmetry-breaking 
state. However, it turns out that this candidate bulk solution is inevitably riddled with 
IR divergences of its flux densities and action, whose physical meaning and implications 
for models of string cosmology call for further investigation. 
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1 Introduction 



Metastable supersymmetry-breaking is an attractive mechanism from a phenomenological 
point of view [1]. Furthermore, theories for which a metastable supersymmetry breaking 
state can be reahzed — such as = 1, SU{Nc) SQCD in the free magnetic phase with mas- 
sive flavours — are relatively simple and generic enough, unlike the comparatively baroque 
ingredients involved in other approaches to dynamical SUSY-breaking (see for instance [21 Ej 
for a review). 

Attempts have been made to embed the proposal of Intriligator, Seiberg and Shih into 
string theory (see for instance |4l |5]), via brane engineering of the electric and magnetic 
phases |6]. Nevertheless, in view of the obstruction that seems to arise upon turning on the 
string coupling Qs ^ or the alternative view [8] that involves string tachyons corrections 
to argue that the brane configuration still describes the ISS state, it is of interest to try 
and find an alternative stringy embedding and search for would-be supergravity duals to 
metastable supersymmetry-breaking states. 

One would (i) start from a well-studied BPS solution of IIB, IIA or eleven-dimensional su- 
pergravity, then (ii) add some supersymmetry-breaking ingredients — typically anti-branes 
at the bottom of a warped throat — so as to lift to a de Sitter local minimum of the poten- 
tial [9] and obtain by the same token a stringy realization of the inflaton potential [10] . The 
next step (iii) is to consider the backreaction of such anti-branes on their background. If this 
can be achieved without any serious singularity or instability, that backreaction procedure 
would then yield a supergravity dual to a metastable state that is part of the same theory 
as the vacuum described by the initial, unperturbed supergravity background. 

For instance, a well studied vacuum of type IIB theory is the Klebanov-Strassler solution 
(KS) [11], which has already less than maximal supersymmetry. It is dual to an A/" = 1 
supersymmetric SU {N + M) x SU{N) gauge theory exhibiting interesting features such as 
a cascade of Seiberg dualities [12], confinement and chiral symmetry breaking. 

A realization of metastable SUSY-breaking starting from the KS background has been 
proposed by Kachru, Pearson and Verlinde (KPV) in the probe approximation [13] (that is, 
step (ii) from the aforementioned guideline). Supersymmetry is broken by adding a certain 
amount of anti-branes which are attracted to the bottom of the throat. These authors 
propose a mechanism in which a fraction of the anti-branes can annihilate (via polarization 
and the Myers effect [2]) with the positive brane-charge dissolved in flux, a process which 
is argued to correspond to the decay of the metastable vacuum in the dual field theory 
description [13]. A related proposal has recently appeared in the work of Klebanov and 
Pufu [15], in an 11-dimensional supergravity context. 

A recent program investigating the construction of metastable states beyond the probe 
approximation (corresponding to step (iii) above) has been initiated in [TB | [T7 t [T5 | [T ^ [2U | 1^ . 
The conclusion from that work is that there is an unavoidable singularity in the IR region of 
the backreacted solutions which have been considered so far. 

In order to determine the nature of such singularities, it is worth understanding whether 
their appearance is related to the particular choice of the background or if it is instead a 
general feature of the backreaction procedure. There are indeed backgrounds which share 
enough features with the Klebanov-Strassler background to be candidate setups for arguing 
about the presence of meta-stable vacua. 

In the present paper we will focus our attention on a type IIA setting, for which we will 



2 



examine the backreaction due to anti-D2 branes. The configuration we will start from is a 
non-singular fractional+ordinary D2 branes supergravity solution due to Cvetic, Gibbons, 
Lii and Pope |22] (CGLPfl 

We find that the candidate IIA supergravity dual to metastable SUSY-breakin^ that we 
build is riddled with singularities arising from the linearized deformation of either the RR or 
NSNS field strengths. Of much concern, those are non-finite action singularities. A novelty 
of our work compared to [Ml Elj is that those singularities are not sub-leading compared to 
the kind of singularities that are allowed as a physically sensible ones, that is those stemming 
from the effect of anti-D2 branes smeared on the S*^ at the bottom of the tip. 

Whereas for the backreaction of anti-D3's on the Klebanov-Strassler solution one could 
have expected, with hindsight, a singularity to arise in analogy with the IIA brane engineer- 
ing of four-dimensional gauge theories, a similar argument does not hold for string theory 
constructions of 2-|-l-dimensional gauge theories. 

Indeed, the profile of the NS5-branes featured in those brane engineering constructions is 
generally not rigid but is instead sourced by the stack of Dp branes in-between (see |25] for 
pointers to the literature and much more on the physics of those brane constructions). For 
four-dimensional field theories living on D4-branes between two NS5's, the profile determined 
upon solving a Laplace equation is logarithmically running. This corresponds to the log- 
running of the gauge coupling for asymptotically free theories. 

On the other hand, for three-dimensional field theories living on D3-branes between two 
NS5's, the profile decays as 1/r away from the location of the D3's on the NS5. Such a mode 
does not have the potential ability to enhance small IR fiuctuations into log-running ones, an 
ability to which one might roughly ascribe the singularities encountered in the holographic 
approach to realizing metastable states in string theory, if those singularities are deemed as 
truly pathological. 

So, proceeding in analogy with brane engineering constructions, for 2+1-dimensional IR 
perturbations should be expected not to affect the UV asymptotics of the background. As we 
shall see as an outcome of our linearized deformation analysis, this is not quite the case for 
the candidate supergravity dual to a 2+1-dimensional metastable state. The IR singularities 
we find are affecting the UV behavior, in the sense that they cannot be completely tamed 
without switching off at the same time the force felt by a probe D2-brane in the UV. 

Besides, having their legs in the wrong directions, those IR divergences cannot be iden- 
tified as the remnant signature of an NS5 instanton through which the metastable state is 
been argued to decay in the probe approximation [131 [U] . 

Such singularities cannot be identified either with those characterizing fractional branes on 
Ricci-fiat transverse geometries before the resolution or deformation of those manifolds (so- 
lutions of the Klebanov-Tseytlin [26] type, whose singularities get resolved in the Klebanov- 
Strassler solution). 

The situation is quite puzzling and it might well be that those singularities are an arti- 
fact of having to smear anti-branes in order to make the problem tractable. Some recent 
results [271 [2H] suggest however that a localization procedure is bound to make things worse, 
rather than alleviating them. Furthermore, an analysis at full non-linear level [29] gives 

^The supergravity dual to the deconfined phase of the underlying theory has recently been considered by 
one of us 

^See 123] for a generalization of the ISS model to lower dimension. 

■^We keep all the options open on this issue. Another comment on this appears in section |6l 
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evidence of how a divergent behaviour is still present even when completely localized sources 
are considered. 

It is therefore a reasonable alternative viewpoint consider that such Coulomb-like diver- 
gencies as physically meaningful, and might somehow be used to discriminate among solutions 
of the string theory landscape. 

We propose the following analogy for linking an eventual singular behaviour and its 
physical causes. In QCD, there are free quarks in the linearized approximation. Their 
"backreaction" results in a Coulomb-like singularity. We know that this is an indication that 
quarks are not good approximations at all to finite-energy states from the spectrum of QCD, 
which instead consists of confined, colorless states. 

Now, the singularities we find involve in particular an IR-divergent NSNS flux-density. 
We would like to suggest, following the above-mentioned situation in QCD, that those sin- 
gularities perhaps hint that some of the scenarios that have been proposed in a probe ap- 
proximation to uplift an AdS-vacuum to a metastable de Sitter one using brane sources do 
not engineer acceptable states of the "spectrum" of string theory. But how is the analogy to 
hold with the Coulomb-like singularities of QCD, given that (anti-) branes are not expected 
to source NSNS flux ? Or do they ? 

Indeed, they certainly do not if we stick to the guideline and intuition drawn from the 
supergravity solutions describing such sources in fiat space. Nevertheless, it has been pro- 
posed that they naturally do [HI [30] in perturbation theory around a complicated, warped 
geometry such as the Klebanov-Strassler solution or the CGLP solution we investigate in 
the present paper. 

Besides, the authors of [TSl ED] argue rather convincingly that the IR-divergencies that 
seem to affect, at linearized order in the SUSY-breaking parameter, the backreaction by 
antibranes of an underlying warped background should disappear at full non-linear order. 
The claim goes as follows: the singularities in the flux densities are naturally sourced at 
linear-order of perturbation theory by the acceptable 1/r behavior of the deformation of 
the warp factor due to, say, the bunch of anti-D3's smeared on the of the Klebanov- 
Strassler geometry. It is then advocated that the 1/r contribution of first-order perturbation 
theor}0, summed up with the l/r^ contribution at second-order, and so on with all the 
other contributions, are nothing but terms in the expansion of the inverse warp factor of 
the backreacted background modified by the presence of antibranes. The whole sum of the 
individual contributions at each order is then claimed to be a perfectly regular quantity, h~^, 
the inverse of the backreacted warp factor. As a result, it is argued that singularities in the 
fluxes are an artifact of perturbation theory that should wash out at full non-linear order. 

It is currently a daunting task to check if this possibility is indeed realized. However, we 
link on this possibility in section [6], where we argue how the issue of those singularities in the 
smearing approximation could maybe be settled by considering 2nd-order expansions for the 
deformation modes of a BPS background, a task which has not been attempted so far. 

But, to come back to the tentative analogy with QCD, it might seem after, in view of the 
above chain of argument involving the full backreaction by antibranes (something we do not 
attempt; we stick at linearized deformation), that this analogy does not hold when applied 
to the backreaction of antibranes at full non-linear order. So what ? 

The reason we maintain that this analogy might possibly hold is that the afore-mentioned 
argument seeking to explain how the singularities in the fluxes should vanish at full non- 

denotes the radial variable in the bulk. 
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linear order is not entirely water-tight. Indeed, a recent paper [22] very convincingly shows 
that in some instance a singularity in the flux density is still present at full non-linear 
order ! 

It is beyond the scope of the present work to offer more credence to vindicate or dispel 
the following possibility but it is very tempting to imagine that the IR singularities we keep 
on flnding upon backreacting the effect of antibranes on some BPS background are similarly 
a hint that some of the constructions which have been proposed has duals to metastable 
SUSY-breaking might instead belong to some "swampland" [32] once the backreaction of 
the SUSY-breaking ingredients is duly taken into account. 

Our results are organized as follows. In section [21 we review the CGLP solution, linearized 
perturbations of which pervade the bulk of the present work. We then recall the basics of the 
method developed by Borokhov and Gubser in section [31 As our analysis makes extensive 
use of a superpotential for the CGLP background, we also outline how the latter is derived, 
as flrst found by Herzog [31]. Next, we explicitly evaluate the force acting on a probe D2- 
brane in section [H For a lighter reading, a key part of our work has been relegated to 
appendix [XI where we expose solutions to one of the system of coupled ODE's governing the 
supersymmetry-breaking modes. As for the remaining modes, their IR asymptotics appear in 
appendix [XI as well. In order to calibrate all that machinery, section [Slpertains to adding BPS 
D2 branes to the CGLP background. Finally, in section [61 we impose IR boundary conditions 
appropriate to the backreaction of anti-D2 branes smeared at the tip of that warped throat 
background. We explain how an IR-divergences in the RR or NSNS flux inevitably comes 
about. 



2 Ansatz for the perturbation 

We start with explaining the Ansatz under consideration in this paper. It respects the 
symmetries of the supersymmetric regular + fractional D2 branes supergravity solution found 
by Cvetic, Gibbons, Lii and Pope. Non-supersymmetric, linearized deformations of the 
CGLP background are part of this Ansatz. 
In Einstein frame, it is given by 

dslo = e-^'^'^ 7]^, dx" dx" + f e^'^""^ \h{r)^ dr^ + e^" (D/i*)' + e^^ dVl^ 

Qs F4 = K{r) d^xAdr + 2m {gi{r) + Cs) J2 A J2 + 2 m {gi{r) + C3) U2 A J2 

+ mg[{r)eijkl/ dr A Dfi^ A j'' , (1) 

iB2 = m [g^ir) U2 + gsir) J2] , = , $ = <l>(r) , 

where dQ4 = e°e" and 

Dfi' = dfi' + e,,fc fi" (2) 

denote line elements on the S^ flber over the S^ base. The coordinates fi\ i = 1, 2, 3 obey the 
constraint = 1. The A^ = A^e"" refer to su{2) Yang-Mills instanton potentials. Their 
fleld strength components 

X = dA' + \ e;., A^ A A"" (3) 
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satisfy the algebra of the unit quaternions, i.e. 



Ji p 



-6'^ 6. 



al3 



^ k '^alS 



(4) 



All in all, this makes the seven-dimensional transverse space a cone over a squashed CP^ [33]. 
The Bianchi identities are satisfied in view of the following definitions and identities: 



J2 =ii'J\ U3 = Dfi' A / 
dU3 = 0. 



dJ2 = Us 



(5) 



2.1 The underlying super potential 

We are next going to reduce to a one-dimensional sigma model the bosonic part of the 
IIA supergravity action. The reason for doing so is as follows: the machinery we rely on 
in order to conveniently obtain the linearized non-supersymmetric deformations around a 
supersymmetric background involves a superpotential. By a superpotential, here we mean 
a convenient book-keeping scalar combination of the fields entering a given supergravity 
Ansatz. By definition, field-space derivatives of the superpotential times the sigma-model 
metric yield the potential that is obtained upon reducing the higher-dimensional supergravity 
Ansatz to the one-dimensional sigma model. 
In Einstein frame, the IIA action reads 



d^^x^/\Y\R 



rf$ A i<d^ + 9se ^ H3 A -kH's 



where 



F^ = F^-CiAHs, Fi = dCs, H3 = dB2, F2 = dC, 
Inserting the above expressions for the fields and metrics ([1]) yields 

£5Vol(Mi,2) Vol (Me) 



(6) 



(7) 



'IIA 



2k 



.2 



dr C . 



where C = T — V and Mi^2, Mg denote the 2+1 dimensional Minkowski space and the level 
surfaces of the seven-dimensional G2-holonomy manifold, respectively. 
The kinetic term is 



T 



-,2u+iv 



h 



30 z'^ + 2u'^ + 12 v'^ + 16 n't/ -2g^ 



-1/2 



m 



-9z+^/2-2u-iv 12 



{g'2'e-'- + 2g','e-'^)-U'' 



2 £6 



(9) 



and, as anticipated, it is clear that both h{r) and K{r) are non-dynamical fields. 

Then, upon eliminating the non-dynamical K through its algebraic equation of motion, 
namely 
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and evaluating the Lagrangian at the corresponding minimum for K, the potential becomes 
V=-2h [e2"+8- - e^-+^^ + 6 6^"+^"] + 2 /i ^ e"^^"* [g^ + ^73]' 

Writing the Lagrangian as 

C = -\Gab (dr/dr) {d<p'>/dr) - V , (12) 
where we denote the set of functions 0"^, a = 1, 7 in the following order 

0" = (m,i;,z,<I>,^i,^2,^3) , (13) 
we find that the following superpotential, initially found by Herzog [31] . 



W = -S [e"+^^ + e'-^'^] + 8 ^ g]/' e'^-t [^^i (^72 + ^73) + 92C, + g, c,] 

correctly accounts for all the terms in the potential ffTTl) . that is to say 

1 dW dW 
V = - G"^-—-— . (14) 
8 90'^ d(t)^ ^ ' 

One can check that the zeroth-order CGLP solutions that we are about to summarize below 
obey the first-order BPS equations derived from this superpotential: 

±G''''^ = 0. (15) 



dr 2 d(f)^ 

This motivates the respective choice of signs in front of the metric part and the flux pieces 
of the superpotential f|T^ . which are otherwise arbitrary. 



2.2 Zeroth— order solution 

The solution of Cvetic, Gibbons, Lii and Pope corresponding to resolved fractional D2-brane 
with transverse seven-dimensional squashed cone over CP^ preserves 1/4 of the original 
supersymmetry, giving rise to a dual A/" = 1 field theory in 2+1 dimensions. 

Let us now gather the expressions for the zeroth-order functions entering this solution, 
around which we will next expand. It might be appropriate to remind that the radial variable 



7 



runs from one to infinitjif]. 



1 7 



-1/2 



2 7;0 -'-2 

4 



1 7 



1 2 

- r , 

2 ' 



9', 



9l 



9l 



h{y) dy, 
f2{y) dy, 
h{y)dy, 



ji — e tti 



/2 = /ie2"^2 

/3 = /ie2^°W3: 



Ml 



Ms 



1 (3r^-l)P(r) 

4r4 - 1) ~ 4^5 (^^4 _ 1)3/2 

1 P(r) 
^ "'^ 1771^^ 1 u/2 ' 





1 



+ 



Pfr) 



2 (r4 - 1) ^ (^4 _ 1) 



3/2 



(16) 



along with 



8z 







2 /"oo 



^/ y [usiy) - U2{y)] ui{y) dy , 



where 



9sHq 



P(r) 



1/4 



(17) 
(18) 



du 



K{-1) - F(arcsin(l/r) 



-1). 



(19) 



From now on, F{(f) \ k) denotes an incomplete elliptic integral of the first kind and K{k) = 
F{n/2 I k). We also will encounter elliptic integrals of the second kind E{(f) \ k). See 
Appendix |B] for a quick reminder. 

The expression for the warp factor Hq abov^ arises from the condition on the trace of 
Einstein's equations 



-(0 1 2 







— m'\G' 
6 



|2 
4 I 5 



(20) 



where we generally define G4 via F4 = K d^x A dr + mG^ and we have used the fact that 
G4 = -k-f Gg for the zeroth-order solution. This can be integrated to 



^10 [e""/' d'x A dH^') = -g;'/' m G° A P! 



I 

2 iMg ! 



(21) 



with Mg the level surface of the G2 holonomy manifold which is part of the CGLP solution 
and ensures its regularity. 



^From now on, superscripts will refer to the perturbation order, while subscripts label different functions; 
quantities which are not labelled by a superscript do not enter the set of perturbed scalars which we introduced 



in UM- 

^It is straightforward to check that the expression PT)) for the warp factor is identical to the one provided 
by Herzog in |31j . i.e. Hq = j^^- y [2 1x3(7/) — 3] ui{y) dy , taking into account different conventions. 
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Note that the UV behavior of the warp factor is as expected, namely 

ifo(r) = ^ + 0(l/r^) , r^oo. (22) 

The IR asymptotics of H is 

7 TO^ 

H,{r) = Ho{l) - - — (r - 1) + 0(r - l)^ r ^ 1 , (23) 

which estabhshes that r = 1 is a coordinate singularity and that the metric is actually regular 
there. Indeed, notice that in the IR the unperturbed ten-dimensional metric takes the form 



dsl, = Ho{l)-'^'ds'+Ho{r 



3/8 



(24) 



_4(r - 1) ' /V r / 2 

The coordinate singularity at r = 1 can be eliminated by shifting gears to 

T = Vf^. (25) 

The space transverse to the D2 branes therefore approaches x S"^ in the far IR. 

As for the constants C2 and C3 appearing in our Ansatz ([T]), the background only specifies 
their difference 

C2 - C3 = ^ . (26) 

Amazingly, it turns out that g^, g2 and g'^ can written in terms of the functions /i and 
/2 or /s appearing in ( fT6l) : 



/i, ^73°=(l + 8e^""+^^")/i. (27) 



^7? = ^e-^"°+^^'V2-C2 or equivalents ^0 = ^e^"" /s - C3 . (28) 
Note that those modes are well-behaved in the IR and their series expansions go as 

^7? + C2 = ^-^(r-l) + 0(r-l)^ g', = l{r-lf' + 0{r-lf\ 

g's = -lir- If" - 4 - 1)'/' + 0{r - Ifl^ , (29) 

which might not be immediately obvious from their defining formulae. Actually, the radial 
derivatives of g^ and g2 also make their way in ^^3. Correspondingly, a piece of the NSNS 
flux behaves as dr A J2 in the infrared. This does not signal a pathological behavior and 
instead is just another instance of the unphysical and non-pathological, coordinate singularity 
dr"^ ~ that we have already encountered. 

3 The Borokhov— Gubser method 

The method proposed by Borokhov and Gubser in [31] allows to find non-supersymmetric 
supergravity solutions, starting from a given BPS background. The idea behind the technique 
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is as follows: rather than having to solve n second-order equations for the n fields (pa entering 
a supergravity Ansatz encompassing the background solution, we trade those complicated 
2nd-order differential equations for 2n first-order equations governing those fields (pa and 
their "canonical conjugate variables" 

The simplicity of the method has much to do with the fact that always form a close 
system. The equations for the modes of much physical interest, (pa involve the conjugate 
modes as source terms. Let us quickly review that approach of Borokhov and Gubser. 

We rewrite the Lagrangian by means of the superpotential (IMI) as follows 

1 1^ dw\ fd4>b 1^ dw\ idw 



The gradient flow equations obeyed by the underlying BPS solutioiill read 

Furthermore, the "zero-energy" condition arising from the Grr Einstein equation is a con- 
straint that applies to any solution, BPS or not: 

-iG*^ + VW = 0. (32) 

z dr dr 

The method of Borokhov and Gubser [M] relies on a superpotential to determine pertur- 
bations to a solution of f l3T]) that satisfy the equations of motion but not necessarily f l3T]) 
itself. Let us consider an expansion of the flelds (pa around their supersymmetric value 0°, 

<Pa = <Pl + <Pl{a) + 0{a') (33) 

for some set of parameters a. Let us introduce the following notation 

e = G'^'m - iV,^(0o)0i) , where N,\(f) = 1 A • (34) 

Inserting the expansion fl33l) into the equations of motion derived from the one-dimensional 
Lagrangian, and keeping terms up to the linear order, one obtains 

^+eN,^{(P') = 0, (35) 
'^l-N^\^')(Pl = Ga,{(P')e, (36) 



dr 

while the constraint (l32l) can be written as 



C— = . (37) 
dr 

The functions are a measure of the deviation from the gradient flow equations (1311) . 
Notice that for a supersymmetric deformation all the vanish. The obvious advantage of 
this method is that one can solve separately for the flrst-order subsystem (135|) and then solve 
for fl36p which are again flrst-order. 



^In the case of present interest this BPS solution is the CGLP solution that we have introduced in Section 
2.2. 
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3.1 The first— order equations for the supersymmetry— breaking de- 
formations 



3.1.1 ^ equations 

We present the system (135|) of first-order equations for the fields which are conjugate to 
the linearized deformations (pa of the CGLP background, in terms of a convenient change of 
variables 

L = (6, ei - 6, 6 + 2 ^4, ^4, ^5, ^6, -^6 + ^7) • (38) 

The above combinations were chosen so as to make the corresponding system of equations 
much easier to solve. We actually managed to find fully analytic expressions for the 
conjugate modes. 

The equations are listed in the order in which we have solved them: 



is 



-4 — — he '''' 2 4 
/6 



^3 



(39) 



64/6 



(40) 



^5 



1 , n„0 ^,.0 15z° 

he 2 4 



Are 



(fe + fr) 



g;6g4«o+6.o+<i,o~^ 



(41) 



1/4 

2/6 



^ g-2u0-4i)0-|(10^0+#'') 



2^i/2/6g2."+4."+9."|^ + e^m2(c2 + 9^1% (42) 



-^(102''+'I'°) 



3/4 



e 4 



24e 



2?i°-4i)''+620+|*'' 



C2 + ^??)e6 



12e 



-2«''+620+|$' 



^*°(C3 + ^??)(e6 + Cr) + 6e^^"^;f (^° + ^,3°)f5 



m 9s 2u°-4v'^+t- 

/6 



;C2^?2° + C3^?3° + ^7?(^72+^?3))e3 



(43) 



9^1' 



-h e 



0_/l„0_15 ,0_*" 



-2«"-4?;' 



3/4;6g.' 



0_|_4„0_,_15 ^0_L« 



+ ^3/4^6gi(12«o+8^,o+3020+*0)|^ 



^;6g4«o+6.o+$o^^^ + (7?)e6 + 4/6e4^"+6^°+*°(c3 + ^??)(e6 + fr) 



9sm\c292 + C393 + 9ii92 + ^^s))^ 



(44) 



el 



^3/4^6g««+4i;0+f ^0+^1^ + 3^3/4^6gi(12n0+8t;0+3020+$0)|^ 



24/^e^"°+«^"+*°(c2 + 9'^)ie + Afe'^'^'^'^^\c3 + + fr) 



+ (7.m2(c2(72° + C393 + 9i{92 + ^?3))6 



(45) 
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3.1.2 equations 



As previously done for the ^ equations, we shift the original (j) to a more tractable linear 
combination d), defined aa 



<Pa = (01, 01 - 2 02, 8 01 + 6 03 - 3 04, 8 0i + 16 02 + 30 03 + 04, 05, 06 + 07, 06 " 0?) • 

(47) 

The set of equations explicitly reads 



1 

20 



i>'i = — h e 



2 20 



6 + 26-20e 
46 + 36-20 e 



"°+^''° 0i-2Oe3""+2^'° 



-60e^ 



(48) 
(49) 



^ 10 



46 + 86 + 6-32 6-806 



1/2 



(76 152°/2+<I>0/4 



4m^ 



2 S's ^ ^ 

+ 2 (7]/^ e^"° [4 05 + (c2 + ^??) (8 0i + 8 02 - 03 



(50) 



(51) 



(52) 



^ ^ g-2tt''-4?)°-3^0/2+3$''/4 "^6 gl520/2+$''/4 / ^ ^4„0 ~^ 

2 gs m? I V 

+ 2 (7^" e^"° U 05 + (c2 + ^7?) (8 0i + 8 02 - 03 



e^^° 6 



^^y^m^e^^'^ (4 05 -(C3 + ^??) 03 



^The inverse transformation is 

05, 2 (06 + '^7) , 2 (06 - 07) ) ■ 



(53) 



(46) 
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0; = _ J_ /, e-2«0-4.0_i5.0/2_$0/4 ^i5.0/2+$0/4 ~^ _ 4 _ 5 
+ 80 £^ g«°+4f'' + 152''/2+*'V4 0^ _ gQ £6 g3«0+2?;0+1520/2+*°/4 

+ 40 (4 (^?° + 93) 05 + 2 (2 + C2 + C3) 06 + 2 (C2 - C3) 07 

-(^2°(c2+^?)+^73°(c3 + ^?))04)]. (54) 

4 The force on a probe D2— brane 

In this section we evaluate the force felt by a D2-brane probing a generic linearized defor- 
mation of the CGLP background. At first glance, the expression for that force might seem 
quite involved. Yet, we will show that using the first-order equations of motion most of the 
terms cancel and the final expression is quite simple, involving only a single mode. This is as 
expected from previous work on the linearized perturbations around IIB and 1 1-dimensional 
BPS solutions [Ml EH EI] • 

Let's then expose the analytic expression we found for the force exerted on a D2-brane 
surveying a generic deformation of the supersymmetric CGLP background. 

We choose a static gauge for a brane spanning Minkowski space directions, without any 
gauge field on its world-volume. The DEI Lagrangian reduces to 

Ubi = -V"^' = -T, e-*/^ gf" V-9oo9u922 = -T, e'^/^-^^^/^ gj'/' . (55) 

The only non-zero RR potential is Cmnpi and the part which gives non- vanishing con- 
tribution is given by 

a, = - /C(r) A dx^ A dx^ , = -Kir) , (56) 

9s dr 

with K{r) given in equation f fTOj) . The Wess-Zumino piece of the D2-brane action thus 
reduces to 

Cwz = -V^' = T,^e^^'-''^{Csh,,,, = -Tp-/C(r) . (57) 
6. gs 

We can now compute the force on a probe D2-brane (from now on we fix Tp = 1). At zeroth 
order we have 

p{0)DBI _ -1/2 TT/ ~<!>0/2~15z° _ 4m^ ^-l/2^-$°/2-152°-2M"-4t;° u ^0 , „ ^0 , ^0/^0 , ^OnI 
^ -9s ^0^ J^9s ' e ' h [025^2 +C3fi'3 +fl'llfl'2 +5'3jJ 



£6 

= -K{r) = l^^;i/2e-*°/2-i5.0-2nO-4.o^ ^^^^o + ^^^o + ^o^^o + 
9s ^ 



and, as further confirmation that everything is under control so far, the two contributions 
compensate each other, as they should. 
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As for the first-order contribution to the force, it arises from 

Tp(l)DBI T7{0)DBll^! I -3/4/1,/ 15 , \ _*^_i5^o 



From these expressions, using the first-order equations fl50|) . fl54l) for 03 and 04, as it happens, 
most of the terms at first-order cancel so that the force on a probe D2-brane reduces to 



p(^^^ ^ p(l)DBI ^ p{l)WZ 

3/4 "■e ?3 



2_X3e2^ 
gs (r^ - 1)3/2 



(5^ 



where we have made preemptive use of the analytic solution for the mode ^3, eq. flHHj) . which 
will be derived in the next section. As an aside, the derivative of the Green's function for 
the CGLP background fl6UI) matches the behavior of the force fISS]) (see [T7j for comments on 
this point). Indeed, allowing only for a dependence on the radial variable, the solution to 

□ G = (59) 

evaluated on the CGLP background is 



G{r) = ci + C2 , - F(arcsin(l/r)| - 1) . (60) 

Vvr^'-l / 

5 Prelims: boundary conditions for BPS D2 branes 

In this section, as a matter of exposing our method before we focus on the candidate back- 
reaction by anti-D2 branes, we first derive the boundary conditions which correspond to the 
modes sourced by a stack of branes placed at the tip of the cone. 

Let us then consider a set of N ordinary extremal D2 branes smeared on the S*^ at the 
bottom of the throat. For the CGLP background, we can explicitly evaluate the Maxwell 
charge 

QcGLp{r) = J^^ e*/2 * F4 = ^^^-^vol(M6)bi(^2 + gs) + c^^?^ + c.g^] . (61) 

This quantity exhibits the following zeroth-order IR behavior: 

QCGLP = , (62) 
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as can be seen from 

gU + 9l) + 0,9^2 + c,9l ^ ^(r - 1)^/^ - ^(r - 1)^/^ + O {{r - 1)^^) , (63) 
using equation f l26|) . 

Within the Ansatz we have been considering, a BPS solution describing the addition of 
N ordinary BPS D2 branes smeared on the S"^ in the IR can be found by shifting 92 and 93 
such that the combination 92 + 93 — which is multiphed by 91 in ( 1611) — does not change: 

92 -^92 + ^^ 93 ^93-^- (64) 
This way, the charge is shifted as 

CGLP QcGLP + ^Qd2 , (65) 

with 

Note that the flux through S"^, 

qs^ = W / ^4 = ^'""^A (^1 + C2)vol(^^) , (67) 

stays unchanged under the shifts (!64|) . while the warp factor shifts as 

H,{r) ^ y /^e-^""-^''" [9'M + ^?3°) + C2^72° + C39I + iV] ^Z/ (68) 

and now is endowed with a singularity of the kind 

H{r) ~ 4% ■ (69) 

This is the expected behavior of the harmonic function for Dp branes smeared on an 
within an otherwise ten-dimensional flat space, which indeed behaves as ^tX-t , where p = 2 
and r = 4 for the CGLP background. 



Let us now see in more detail how this BPS solution can be reproduced by the first-order 
perturbation apparatus. First of all, we set to zero all the modes related to supersymmetry- 
breaking, namely we impose that all the constants Xa and Bi ~ X3 fl90|l . which enter upon 
integrating of ^ equations, should vanish. 

Furthermore, the zeroth-order combinations e^"° and e^^° reach constant or zero value 
in the IR ; since we expect that the geometry of the transverse space is not affected by the 
addition of BPS D2 branes we impose the perturbations associated to u and v to vanish as 
well. This fixes 

yiR = y^//? = 0. (70) 
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In addition, non-singularity of 05 and 07 (we recall that they enter the fluxes of our 
Ansatz) is ensured by 

F/« = -^K{-1){7Y^^ + 80F/«). (71) 
The mode 06 is regular, and in view of the first-order contribution to fl6T|) 



32 16 y " 16 " 

(72) 

one should impose that 06(r — ?■ 1) be proportional to the number of BPS D2 branes spread 
over 5"^ at the tip. 



To recap, the above choices of integration constants f l70l) - (17Tl) yield the expected behavior 
for BPS D2 branes added in the supersymmetric CGLP background: 

01 = 0, 02 = 0, 05 = O(r-l), (73) 

= ir/« + o((r - , 0, = + o((r - 1)1/2) 

We recall that 0i 2 denote perturbations of the stretching functions, 03^4 label perturbations 
of the warp factor and dilaton, whilst 05,6,7 are the modes corresponding to the linearized 
perturbations of the NSNS and RR fluxes of this IIA background. 

The integration constant Y"/^ is the only remaining one and is related to the number A^ 
of added BPS D2 branes: indeed, the equations for 0e and 07 reproduce the shift flMl) . The 
warp factor, along with the dilaton, acquires the expected singularity and 

H = e8-« (1 + 8 03) , e* = e*« (1 + 04) = e^^" (1 + 2 03) , (74) 

in accordance with ~ H^f^. 

6 Assessing the anti— D2 brane solution 

The final step and main aim of our analysis is to determine how, within the space of generic 
linearized deformations of the IIA CGLP background, one can account for the backreaction 
due to the addition of anti-D2 branes smeared on the at the tip of the warped throat. 

As the prime physical requirement we should impose that the force felt by a D2 brane 
probing the backreaction due to this stack of anti-D2 branes be non-vanishing. So, we are 
forbidden from turning off the corresponding mode which appears in the expression flSSl) of 
the force, and enters the various expressions for the modes 0^ by means of the shorthand 
combination 

B, = ^X3e-«^o(^) . (75) 
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As our next set of IR boundary conditions, let us recall that the modes 03 and 04 as- 
sociated to the perturbation of the warp factor and the dilaton must exhibit no worse than 
a 1/ y/r — 1 ~ 1/r behavior (cf. equation fl25l) ). Such a behavior is in accordance with the 
Coulomb-like divergence associated to anti-D2 branes smeared over the S'^ at the tip of the 
warped throat. 

Inspecting the IR expansions of the deformation modes 0a, every piece that is more 
singular than the aforementioned 1 / i/r — 1 behavior will be culled by tuning appropriate 
combinations of the X's and the F' s integration constants parametrizing the space of generic 
linearized perturbations of the CGLP background. 

Another, equivalent but slightly less liberal, criterion that we are about to consider fo- 
cuses on allowing or discarding various pieces from the 0a's IR expansions depending on their 
contribution to the energy. More precisely, we consider the kinetic energy ^ and the poten- 
tial energy f|TT]) obtained by reducing the IIA supergravity Ansatz ([1]) to a one-dimensional 
sigma model. 

For instance, the energy associated to the first-order perturbation of the dilaton and warp 
factor is obtained by expanding to second-order the corresponding terms from dH]): 



g2 [u°+(t,i)+A[iP+4,2) 
h 



30 (^°' + 0'3)'-^($°' + 0l)^ 



h 



30 0'3' - \ 0'/ - 2 (01 + 2 02) ($°' 01 + 60 0'3) ] (76) 



The energy associated to the deformation of the warp factor and dilaton exhibits the following 
singular behavior 

/ -|n3/2 / "'/^3,4 



1 

dr J (r — 1)'^/^ ' 

where as a matter of course we neglect less diverging terms. This behavior sets the threshold 
for what we consider an allowable singularity in the energy. 

Note that, as it turns out, for all practical purposes we can neglect contributions of the 
type 0a0fe and 0^ 0b for a ^h: they only contribute to sub-leading divergences. In addition, 
there is no contribution to the energy that is first-order in the SUSY-breaking parameters, 
since we are expanding around a saddle point. 

Another remark is in order. We have considered linearized deformation for the fields 
entering the supergravity Ansatz (JT]), namely we have expanded as 

0, = 0° + 0^(X,r), (77) 

with Xi and Yi being implicitly the small supersymmetry-breaking expansion parameters. 
On the other hand, we are considering quadratic contributions of the 0^'s to the energy. 

The reason why we do not stop at first-order contributions to the energy from those 
deformation modes is that we have expanded around a saddle point. Had we gone as far as 
computing 2nd order expansions of the deformation modes, namely 

0a = 0° + 0i(X, Y) + 02 (X, y, Z, W) , (78) 
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which is an achievable if strenuous task, it might well happen that the singularities we are 
about to expose might cancel against truly second order contributions to the energy. By this 
we mean contributions of the type (/'a^b, in addition to those of the form (0^)^ 0° that we 
presently consider. 

Everything is now in place to show that the candidate IIA supergravity dual to metastable 
supersymmetry-breaking that would be obtained out of backreacting D2S spread over the 
5*^ in the far IR of the CGLP background comes with an irretrievable IR singularity. Indeed, 
we are going to show that it is not possible to simultaneously satisfy the two previously 
mentioned physical requirements. 

In point of fact, there is a singularity associated to the NSNS and RR fluxes that is 
worse than the ones we allow, namely those that are physical and should be kept based 
on their identification with the effect of adding anti-D2 branes to uplift the AdS minimum 
of the potential. There is only one way of getting rid of that "unphysical" singularity: it 
entails setting to zero the single mode entering the force felt by a brane probing the non- 
supersymmetric backreaction by D2S. So, our two sensible IR boundary conditions are 
incompatible. 

Ensuring that there is a force exerted on a probe D2-brane by the anti-D2's at the tip 
results in a ■^"[yr ~ singular contribution to the energy, stemming from the NSNS or the 
RR field strength. Such a singularity is worse than the ones it is sensible to a priori allow, 
namely j^/a singularities or milder ones, associated to the smeared D2^- 

Let us see how this comes about with full details. First of all, note that the potentially 
most divergent deformation modes is (p-j: its IR series expansion f ll27p displays and 

pieces. That mode, 07, contributes only to the deformation of the NSNS 3-form field 
strength 

i6H3 = m [(06 + 07) Us + 0; drAU2 + 0^ dr A J2] . (79) 

In view of ([9]) and (fTTI) . the leading contribution to the energy from the deformation of the 
NSNS 3-form is 



2 ^2u"+iv" 



m e 



h 



6'2e-4«o + 2 0'2e-^^'° 



-2^/^e-«^"[0e + 0rf . (§0) 



There is another potential contribution from ffTTl) which involves 06 and 07. It is easily seen 
that it is sub-leading. Now, what is the IR singular behavior of flHOj) ? We focus on the most 
singular piece of 07 ~ ^^"^ derivative. It entails the following singular behavior 



2 

^ p-8 20{r) 
"£6 



h{r) \dr (r — 1)/ \(r — 1 



1 



- 1)5/2 



According to our physical criterion pertaining to the energy, we should then discard the 
most IR-divergent piece of 07, see (11271) . This is achieved by imposing 

^5 = ^^^ [3 (l7 + /i(-l)2) 5i + 56A'(-l) (3X6-2X7)] , (82) 
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Xi = -J— [6048 Bi + 1032192 X4 + 215040 F/^ - 2580480 YI^ + 215040 Ei-l) F/^ 
86016 L ^ ^ V / ^ 

+ 235200 K{-1) Xq - 139360 A'(-l) Xj + 133120 i^(-l) F/^ 

- 21504 A^(-l) Yi^ - 245760 K{-1) F/^ + 8364 K{-lf Bi 

- 27216 A^(-1)=^X6 + 11304 A:(-1)3X7 - 1809 A:(-l)^5ij , (83) 

where (I82p has been apphed to obtain (1831) out of the combination of X's and F's from the 
P^'^t of 07's IR expansion. 
We now turn our attention to getting rid of the singularities stemming from the RR flux 
and 05. First of all, note that the condition f l82|) washes out, at no extra cost, the leading 

log(r-l) 



/r-1 



part of 05 's IR asymptotics. 



Still, one should enforce that the X-^ part of 05's IR expansion be wiped out by ap- 
propriately tuning some of the X's and F's. Indeed, if kept unchecked, that divergent piece 
would yield a singularity in the energy arising from the RR flux: 

2 -8z0-9<)i3+</'4/2 

_ 4 !^ g~8.o-903+*4/2+2.o+2</,, ^ J2 [gl + C2 + 05)' e"^''"-^'^^ + [gl + C3 + 05)' e-^"°-^^ 

1 



^.5/2 



^4) 



which is beyond the energy threshold fl77|) and should be culled. To get rid of that singular 
piece from 05, one must exact 

- 32 A:(-1) (6384 Xe - 3711 X-j + 4160 Y^^ - 672 Yj^ - 7680 F/^) 
6 A:(-1)^ (1795 Bi - 3976 X5) + 152 A:(-1)^ (336 - 179 X^) 

+ 2235 A'(-l)^ Bi - 42024 Bi + 1344 [64 Xi - 768 X4 - 23 X5 

- 160 (F/^ + E{-i) F/^ - 12 F/^) ] = . (85) 



We have finally reached the punchline of our analysis: taking into account the conditions ( l82l) - 
(|83!) that did arise from ensuring that no "unphysical" singularity pops out of the NSNS flux, 
it turns out that (185!) yields 

11340 (4 - A:(-1)^) 5i = , (86) 

in blatant opposition to the physical requirement that a D2-brane probing the non-supersym- 
metric deformation of the CGLP background experiences a non-vanishing force ! 



We have therefore come to the conclusion that a careful analysis of the backreaction of 
anti-D2 branes on the CGLP background inevitably results in an IR singularity. By focusing 
on two particular flux elements for which the energy contribution can be easily calculated, 
we have shown that it is not possible to avoid a singular behavior provided we want to keep 
the Bi mode entering the expression for the force (!58l) to be non-vanishing. 

One has to face that at least one of the perturbed NSNS or RR fluxes contributes to a di- 
vergent energy density and to a divergent action as well (given that the factor ^/gio ^ ^/r — 1 
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appearing in the ten-dimensional action (jS]) is not enough to make the action finite in the 
IR), much as is the case in [21]. The key difference from fIT\ lies in the fact that in our case 
the singular behavior is not at all sub-leading. 

The above type IIA analysis completes the program of investigating the would-be backre- 
acted supergravity duals to metastable supersymmetry-breaking vacua, which was originally 
started in a type IIB setting [16], and next considered in [21] in an 11-dimensional context. 
It would be of much interest to consider other backgrounds and/or, as explained at the begin- 
ning of this Section, to go to higher-order in the perturbations around those BPS solutions. 
It might be that an absence of the nasty singularities we have kept on encountering so far 
could be used in order to discriminate among solutions of the landscape string theory vacua. 
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A The space of linearized deformations 

A.l The system 

A. 1.1 Structure of the solutions 

We start with some remarks on our approach to determining solutions to the set of equa- 
tions (l39l)-(H5l). In the case at hand, we were able to find fully analytic expression^ for all 
the C,a- We present the comments in the order in which the corresponding equations have to 
be solved. 

The first equation we have to solve is the one for ,^3. Upon recognizing the algebraic 
expression for K[r) ffTOj) and keeping in mind that for a BPS solution the following identity 
holds 

^(0)W = -tA^, (87) 



it can be expressed as 

The solution obviously is: 



77/ 

-Ho 



|3(r) = X3i7o(r)e-«^"(^) . (89) 



^It is important to have a solution expressed in terms of the least possible number of nested integrals. 
As happens in previous similar work (TBI HH [T^ , it is usually not possible to find a fully integrated solution 
and one then has to be content with series expansion ; if the number of nested integrals is important, that 
quickly becomes burdensome. When counting of nested integrals we do not take into account the one which 
enters the definition of the elliptic functions. 
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Let us introduce at this stage the constant Bi, which we find convenient to use in order to 
avoid extra clutter 



5i = ^X3e-«^°« . (90) 

The next step is to exphcitly perform the integration entering the expression for the 
CGLP background warp factor Ho{r), which we rewrite here as 



2 /-oo 



Ho{r) = ^ / r My) - My)] My) dy (9i) 



The integrand has the following structura^i 

a2J^{rf + aiJ^{r) + ao, (94) 

with ttj some functions of r which do not involve J-". We simply apply integration by parts 
(in the following we drop the radial dependence for ease of notation): 

j a^T'^ + axT + = Ag-F^ + J (ai - 27' A2) + j ao 

= A27^ + As7 + J (ao-rAs) 
= + A37 + , (95) 

where the labels introduced above denote the following: 

jr' = Ai?(arcsin(l/?/)| - 1) ^ 



dy y^M^ 
"3 = "1 - 2 ^"'742 , 
= Qq — J-'A^ , 

Ai = j ai. (96) 

Once we have a primitive we have just to evaluate it at the two extrema of integration to get 
an analytic expression for Hq, and therefore for ,^3. 
The equations for ^7 is: 

^^-"647^^' - "4 (r^ - 1)3/2 (97) 

which can be directly integrated. 



^"Wc adopt here, and for the remainder of the paper, the following calligraphic notation for the incomplete 
elliptic integral of the first kind F: 

F{r) = F(arcsin(l/r), -1)) (92) 

and, similarly, later on we will refer to 

£{r) = £;(arcsin(l/r),-l)) (93) 
as the incomplete elliptic integral of the second kind E. Cf. also Appendix [BI 
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The functions ^5 and are coupled into a subsystem of ODE's, which we can rewrite as 
a = -2h {26'-°-'^" + e-'-") - 2/.e-^"" - f /i (98) 
l; = -/^e5-^^e-2""+'^V2el (99) 

In order to obtain a solution, we first have to solve for the homogeneous system ; we arrange 
the two basis vectors of the space of homogeneous solutions in the so-called fundamental 
matrix 



or —1 X'/' 

The solution to the inhomogeneous system is then expressed as 

( Iw ) =^-W^- + ^W / '^^M-'9l{y)dy 

where X55 = {X^,X^ are integration constants, and = {g5,gf) is a book-keeping for the 
non- homogeneous terms entering equations (pS]) and fl^I?]) . 

The equation for ^4 is entirely non-homogeneous and depends on ^5 and ^q. We can 
rewrite it as follows: 

C = Ih h h -\U2 + h) - \ Ml - e-\2us - 3)u^ . (101) 

which we managed to integrate. 

Finally, the functions ^1 and ^2 are entangled into the following system of first-order 
differential equations: 

= he-""' + /ie"°-2-° 6 - 2(/2 - me + 2M7 - ^r(2w3 - 3)u, , (102) 

o 

= h e""" + 3/1 e""-^''" 6 - 2(3/2 - /s) + 2/3 + ^^(^^3 - 3)ni . (103) 
whose fundamental matrix H,2 reads 



2r" 



E 



(1 - rVf43I(g(^) _ jr(y)) 



-2r^ (^(r) - J^(r)) 



(104) 



Analytic expressions for ,^1^2 are listed in the next subsection, along with solutions for their 
siblings. 

A. 1.2 Fully analytic expressions for the modes 

Here, we collate analytic solutions we derived for the ^ syste 



^^Wc made sure that those solutions are exphctly real, which straightforward if gruelling successive inte- 
gration by parts do not immediately yield. 
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112r5(r4 - 1)3/2 

/ 189."- 258.» + .' + 48 ^ + ' \ 

\ 1792r4(r4-l) v i v ; 6 ^ ' 2Q%%r^{r^ - Ifl^ ) 

{ 69ri2 - 114r*^ + eir^ - 24 SlSr^^ - 390rS - 53r^ + 120 

+ ^(r) ^ - 5i 896^3(^4 _ 1)3/2 (-1) 3584^4(^4 _ i) 

^ 4 ^ 63ri2 - 78r8 + 31r^ - 8 Qr^^ - ISr^ - 7r^ + 8 
+ A2(r — 1) — Ae „. 4. 4 TT -^7 



64r4(r4 - 1) 9Qr\r^ - 1) 



+ (24X5 + J^(-1)(24X6 - I6X7 - 35iJi (-1))) 



384r5(r4 _ 1)3/2 



^ Slr^ - 75r^ + 16 „ , 2>lhr^'^ - SlGr^ + 229r^ - 60 , 4 



512r4(r4-l) ' r ' ' ^16r4(r4-l) 



r- + 1 _ 33rS - 35r^ + 4 63ri2 - 78rS + 23r^ - 4 



K,4 



- X5i^(-1)-— - - Xe— — + X6ir(-1)- 



16r5(r4 - 1)3/2 64r3v/;^4^ ^ ^ 64r4(r4 - 1) 

I ^ 9r«-llr^ + 4 ^ 9r^2 _ ^g^s ^ ^4 ^ 4 
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i2 = Hrf[B ^ 



112r5(r4 - 1)3/2 

189ri6 _ 438ri2 + 2Alr^ + 52r^ - 16 {45BiK{-l) - IGSXe + 112X7)(r4 - 3) 
896r4(r4 - 1)2 2688r5(r4 - 1)3/2 

, , , 69ri2 - 132rS + 2br* + 20 , , SlSr^^ - TSOr^^ ^ 427^8 ^ 75^4 ^ 44 

-^W -^1 ..o^...^4 5iA(-l)- 



443^3(^4 _ 1)3/2 V I 1792r4(r4 - 1) 

4 63r^2 _ 87^8 ^ 4Q^4 _ ^3 Qr^^ - gr^ - IGr^ + 12 
~l~ ^2 2r — 7^ — ^ — X'j — — ^ 

4 q 

+ {K{-l){?,BlK{-l)-2AXQ+lQX^)-2AXr 



2 



5y 



384^5(^4 _ 1)3/2 



^ 51r*^ - SOr^ - 32 ^ , SlSr^^ - SGlr^ + 40r4 + 134 , . , , 

89frV-l) ^ 1792.^^ - 1)3/^ + ^' - ^('■) 



256r4(r4-l)2 n6r4(r4 - 1) ^ ^ ^ 16r5(r4 - 1)3/2 

33rS - 38r^ + 6 , , 63ri2 - Slr^ + 32r4 - 6 

' ^2r^V^^^ ' ^ ' 32r4(r4-l) 

9^8 _ 14^4 ^ g g^l2 _ 9^8 _ 8^4 ^ g 

+ X7 +XjK{-l) — ^ , 106 



e3(r) = X3e-«^°(^)i/o(r) , (107) 
where 

_ / 3 1 



2£6 ^ ' \32 8r4(r4- 1)2 

/ 3r^ + 3r^ - 4 K{-1) 

~ W^^' Vl6r3(r4- 1)3/2 + 16 r4(r4 - 1)2 

/ 3r^ - 4 3^8^3^4,4 i^(-l)^ . 

^ 2^ V32r2(r4- 1) + 16r3(r4 - 1)3/2 ) " 8r4(r4- 1)2 ' ' ^ 
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3, 3B,{3r^-l] 



448r5(r4 - 1)3/2^ 
5i(lllr^2 _ 222r8 + 99r^ - 16) 
3584r4(r4 - 1)2 



3584^5(^4 _ 1)3/2 



;i68X6 - 112X7 - 455iA^(-l)) 



fii(15r8 - 12r^ + 10) BiK(-l)(201r^^-402r^ + 4:5r^ + U) 
•r{r) I . tt^ta^ . — \- 



896r3(r4 - 1)3/2 



7168r4(r4 - 1)^ 



- 1 9r^ - 9r^ + 4 

512r^(r^ - 1)3/2 ("24^5 + K{-1) {3B,K{-1) - 24Xe + I6X7)) + ^28^4(^4 - 1) ^^^^ " ^X^; 

5i(51r^-32) 5iir(-l)(201r8-231r4 + 134) 5iJs:(-l)2(9r4 - 5) 3Jr(-l)(3r4 - 1) 



3584^.2(^4 



7168r3(r4 - 1)3/2 



512r4('r4 



- 2 , 
+ ; 3X6 - 2X7 



3X5 + i^(-l)(3X6-2X7) ^ ^ 



64r4(r4 - 1) 



64r5(r4 - 1)3/2 " 5 
(109) 



fii(15r^ 



Bi{l-3r^' 
5iK(-l)(3r^ 



1) (3r^-l)(3X6-2X7) 35i(5rS-5r'^-2) 



8r4(r4 - 1) 3r4(r4 - 1) 

21r^+10) 3BiK{-l) (3r^ - 1) 



28r3Vr4 - 1 



28r2(r4 - 1) 



1) 



X5 + 6r - 1 Xe 



3X5 — 2X7 
3rVr^ - 1 



^6 = Hr) 



Bi 



, _ , -5i(15r8 + 3r4-4) Siir(-i: 



112('r4 - 1) 



3fiir(5r^ + 4) 5ii^(-l) 



112v^ 

f7(r) = X7 + 



1) 



X, 



r-\/r4 — 1 



3Xfi - 2X' 



1 Xr — — X7 , 

2 / ^ 3 ' 



(110) 
fill) 



A. 2 The 0a equations 

A. 2.1 The structure of their solutions 

The system does not admit a fully analytic solution. We are thus forced to consider either 
numerical work as in [20] , or to series expansions, that latter option meeting our the needs of 
the present paper. Here, we present the equations, and show that we can find solutions up to 
three nested integrals. This might prove helpful to future work. Once more, the presentation 
follows the order in which the equations have to be solved. 

We report our final results back again to the (p basis, as it is convenient to impose boundary 
conditions singling out the effect of anti-D2 branes in these variables. 
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The functions 0i and 02 are coupled and the system is 

1 



0'^ = -/i e"" 01 - 3/i e 
The corresponding fundamental matrix is 



20 



(4^1 + 36 



:il2) 
:il3) 



3-r* 



r^Vr''— 1 



^ + ;37;fe(^W--F(r)) 



3-ri 



-r 



(114) 



A formal solution is thus 



0i(r) 
Mr) 



(115) 



where F12 = (Yi, K,) are integration constants, and gffa = ((^f , (7^) encodes the non-homogeneous 
terms in the couple of equations (I112l) - flll3p above. Some of the integrals can be explicitly 
done but sadly there are some terms for which we were unable to find a primitive. We thus 
have a semi-analytic solution, that is up to an implicit integral. 

We can use the following relation arising from the equation for 0i, 



-he "01-/16 



h 

20^ 



-2u0-4i;0 



(116) 



in order to simplify the equation for 03, which will then take the form 



and has the following solution 

03(r) = 80i(r) + ; 



10 



6 



(^4 _ 1)3/2 



dy 



6-326 



256 



:il7) 



6 



(^4 _ 1)3/2 



dy + Ys 



:il8) 



which is again implicitly defined in terms of a single integral. 

As for the modes 05 and 06, they are coupled and the relevant sub-system is 



>1 - <P3 
£6 



^5 ^ 



h HqC 



2u"-Av'- 



6 + 



-2m'- 



2m2 



4 



h + 802 

whose fundamental matrix is 



/3, 



1 



r4(r4-l) 



' ^ 2 + 3r^) + j^^(r) 



21 



2{6r^-6r*-l) . 2(l-3r"^) 
21r-3VPl-T 21r4(r4-l) 
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A formal solution will have the same structure as (11151) . Recall that gf^ features quantities 
defined in terms of one implicit integral coming from 0i, (f)2 and 03. Consequently, the 
expressions we get are defined in terms of two nested integrals. 

The equation for (pj can be cast into the form 

4>7 = 0'6 - Hoe-'''' + ^/303 - 4/i°e-2"V5 , 
where /s features in equation f|T6l) . and its solution is given by 

Among the summands which appear under integral sign, the first contains no further 
integral whereas the second integrand is itself defined implicitly and so counts as two nested 
integrals. The last summand involves three nested integrals (one explicit here and two coming 
from 05). A simple integration by parts can reduce that number by one, which results in an 
expression for 07 that contains at most two nested integrals. We obtain 

Mr) = Mr) - —^hHoe-'-° ^r) + \hMr) 

We can now use the 0i, 02 system to simplify the equation for 04 that is obtained 
from fl36l). which can be recast to 



61 = -Fo-^i7^04 + 160; - 80', + -^he-'-'-'^'^l, - ^-^h H^'e~'-"~'^" 



- ^^-"°^o-V206 + l^e-'^^'-'^'kH,' (0e - 07) . (119) 

The homogeneous solution to this equation is (pA^hom = . Labelling by gf the non- 
homogeneous piece of flll9p . a general solution is given by 

04(r) = H^\r)Y, + H^\r) Hoiy)gt{y)dy . (120) 
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A. 2. 2 IR asymptotics of the 0^ modes 



We collect here the IR expansion of the (pa fieldo- We write explicitly only the most divergent 
and constant terms, since higher order terms of the IR expansions (we recall here that the 
far infrared corresponds to the limit r — )■ 1 in our conventions) do not provide any constraint 
on the space of solutions. We also impose throughout the zero energy condition flTTj) which 
requires that X^^ = 0. 



4480 



3Bi (^34 + 65A'(-l) 



1792X1 + 336X5 - 112/C(-1)(^3X6 - 2X7) j + 0(^(r - ly/^^ 
3Bi (a1 + lOOK(-l)^) + 2688X1 + 924X5 



fl211 



13440v/f^ 



- 308i^(-l)(^3X6-2X7 



1 



+ 1)1/2 



(122) 



15482880 {K{-iy - 4) 
- 56(K(-1)(2X7 - 3X3) + 3X5 



4801og(r - 1) {K{-lf - 4) (3B1 {K{-iy + 17) 



42i^(-l)=^(21067fii - 49152X4 + 17384X5 



+ 87369^1 i^(-l)^ - 374856^1 - 32i^(-l) (^189168X6 - 1201 17X7 

+ 32(5210F/^ - 33(7K/^ + 80^6^^))) + 40/s:(-l)^ (36624X6 - 22535X7 



+ 1344(^4160Xi - 19200X4 + 311X5 + 160(7r/^ + 7^2 £^(-1) + 132^/" - 144^7^ 

+ -^(35ii^(-l)-8X7)log(r-l) 
256 



2Y, 



IR 



3(K(-l)2-4) 96 



l(48F/« + F/«)+0(( 



a/2 



(123) 



^They can are easily obtained from the (f>a modes via the inverse transformation p6)) . 
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480 log(r - 1) {K{-iy - 4) (sBi (/i (-1)^ + 17) 



6/i (-l)2(92035i - 56(92160X4 + 678IX5)) 



7741440 {K{-iy - 4) 

- 56(K(-1)(2X7 - 3Xq) + 3X5) 

- 301355iA'(-l)^ - 2254920^1 + 32 (^if(-l) (488208X6 - 331467^7 

+ 32(-5210F/^ + 231Y3^^ + 2640F/^)) + 42(4160Xi + 79104^4 + 4919^5 

+ 160(7F/^ + 7F/^^(-l) + 132Y5^^ - 1441^7^^))) + 8A'(-1)^(338909X7 - 494256X6) 
1 



+ — (3fiiA'(-l) - 8X7) log(r - 1) - 



128 



4^4^^ 



2 16 



o((r- 1)1/2 



(124) 



60 log(r - 1) {K{-iy - 4) (sBi {K{-lf + 17) 



5160960 y/F^ 

56(A:(-1)(2X7 - 3X6) + 3X5) j + 6A:(-1)2(1795Si - 3976X5) + 2235SiA:(-1)^ 

42024^1 - 32A:(-1) (6384X6 - 37IIX7 + 4160F/^ - 672K/^ - 7680^6^^) 

152A'(- 1)^(336X6 - 179X7) + 1344 (64Xi - 768X4 - 23X5 - 160(lf^ + Y^^E{-1) 



12F/«)) 
1 



-(35ii^(-l) - 8X7) + 0(r- 1)1/2 



2048 



6A'(-l)2(36599Ei + 3O856X5) - 5115EiA'(-l)^ + 140376Ei 



(125) 



20643840 

1344(l262Xi - 2304X4 + 185X5 + 160(5F/^ + 12^5^^ - 48Y7^^ + 5r/^A;(-l))) 
32A:(-1) (28560X6 - 18495X7 + 44480F/^ - 672Y3^^ - 7680^6^^) 



8A:(-1)^(16841X7 - 26544X6) 



+ 0((r-l)) 



(126) 
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6K{-lf {5656X5 - 2955i) - 8fs:(- 1)^(7644X6 - 4241X7 



5160960(r - 1) 

- 24155iA'(-l)^ + 32K(-1) (7644X6 - 4551X7 + 32(130r/^ - 2iy3^^ - 240^/^)) 
+ 8904^1 - 1344 (64Xi - 768X4 + 7X5 - 160(F/^ - 12Y'/^ + F/^E(-1))) 



log(r- 1) K{-lf -A 
(r- 1) 86016 

log(r — 1) 



35i(17 + K{-lY) - 56(3X5 - J'S:(-1)(3X6 - 2X7)) 



860160 



5i(3468 + 4485i^(-l)2 - 15K{-1)^) - 56(768Xi + 204X5 



- 68A'(-1)(3X6 - 2X7) - 15A'(-1)2X5 + 5A'(-1)3(3X6 - 2X7)) 



+ 



20643840 



32A:(-1) (28650X6 - 18495X7 + 32(1390F/^ - 21K/^ - 240Y"6^^)) 



+ 6A:(-1)2(365995i + 3O856X5) - 8A'(-1)^ (26544X6 - I684IX7) - 5115A'(-l)^5i 
+ 140376^1 - 1344(1216X1 - 2304X4 + I8IX5 



+ i60(5y/^ + i2y5^^ + 48y7 + 5y/^E(-i))) 



B Elliptic functions 



+ C((r-1) 



1/2 



(127) 



For the reader's convenience, we list in this section the definitions of elhptic functions of which 
we make frequent use in the bulk of the text and especially in Appendix A. The incomplete 
elliptic integral of the first kind is defined as 



F(0, m) 



d^ 



'0 a/i — sin^ 6 
while the complete elliptic integral of the first kind reads 



:i28) 



K{m) = F (— , m) . 



(129) 



Analogously, E{(f),m), the incomplete elliptic integral of the second kind, has the following 
expression 

E{(l),m)= I Vl -m2 sin^^rf^, (130) 
Jo 

and is related to the complete elliptic integral of the second kind as 



E{m) = E ( — ,m) 
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